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^' Abstract 

^^ ■ We study a diffusion model of phase field type, which consists of a system of 

^O , two partial differential equations involving as variables the thermal displacement, 

I that is basically the time integration of temperature, and the order parameter. 

Our analysis covers the case of a non-smooth (maximal monotone) graph along 

with a smooth anti-monotone function in the phase equation. Thus, the system 

turns out a generalization of the well-known Caginalp phase field model for phase 

r^ • transitions when including a diffusive term for the thermal displacement in the 

Cd ■ balance equation. Systems of this kind have been extensively studied by Miranville 

^ I and Quintanilla. We prove existence and uniqueness of a weak solution to the 

initial-boundary value problem, as well as various regularity results ensuring that 

CN ' the solution is strong and with bounded components. Then we investigate the 

^1 , asymptotic behaviour of the solutions as the coefficient of the diffusive term for the 

(y-s I thermal displacement tends to and prove convergence to the Caginalp phase field 

0^ ■ system as well as error estimates for the difference of the solutions. 

cn ; 
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^ . 1 Introduction 

This paper is concerned with the initial and boundary value problem: 

Wtt-aAwt- pAw + ut = f infix(0,T) (1.1) 

ut- Au + -f{u)+g{u)3wt infix(0,T) (1.2) 

dnW = dnu = on r X (0, T) (1.3) 

«;(■, 0) = m;o, Wt{-,Q)=VQ, u{-,0)=uo inQ (1.4) 
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2 Solvability and asymptotic analysis of a phase field system 

where ^2 C M^ is a bounded domain with smooth boundary F, T > represents some 
finite time, and dn denotes the outward normal derivative on F. Moreover, a and f3 are 
two positive parameters, 7 : R — )■ 2^^ is a maximal monotone graph (one can see [21 in 
particular pp. 43-45] or [1]), (7 : R — )• M is a Lipschitz-continuous function, / is a given 
source term in equation (11.11) and wq, vq, uq stand for initial data. The inclusion (in place 
of the equality) in (11.21) is due to the presence of the possibly multivalued graph 7. 

Equations (ll.ip - (ll.2l) yield a system of phase field type. Such systems have been 
introduced (cf. [3]) in order to include phase dissipation effects in the dynamics of moving 
interfaces arising in thermally induced phase transitions. In our case, we move from the 
following expression for the total free energy 

^{9, u)= I (-U^ -9u + </)(n) + G{u) + ^|Vnpj (1.5) 

where the variables 6 and u denote the (relative) temperature and order parameter, respec- 
tively. Let us notice from the beginning that our w represents the thermal displacement 
variable, related to 9 by 

w{-,t) = wo + {l*9){-,t)=Wo+ l9{-,s)ds, te[0,T]. (1.6) 

Jo 

In (11.51) . (f) : [0, +00] — )■ M is the convex and lower semicontinuous function such that 
0(0) = = min0 and its sub differential d(j) coincides with 7, while G stands for a 
smooth, in general concave, function such that G' = g. A typical example for and G is 
the double obstacle case 

0(n) = /[_i,+y(n) = r '!!''! -|, G{u) = l-u' (1.7) 

1+00 it \u\ > 1 

so that the two wells of the sum 0(m) +G{u) are located in —1 and +1, and one of the two 
is preferred as minimum of the potential in (11.51) according to whether the temperature 
9 is negative or positive. Indeed, note the presence of the term —9u besides (f){u) + G{u) 
in the expression of \&. 

The example given in (11.71) is inspired by the systematic approach of Michel Fremond 
to non-smooth thermomechanics: we refer to the monography [7] which also deals with 
the phase change models. In the case of (ll.7p the sub differential of the indicator function 
of the interval [—1,-1-1] reads 



^ G 9/[_i^_i_i](n) if and only if ^ 



Let us point out that, with a different terminology motivated by earlier studies on the 
Stefan problem [6], some authors (cf. [7j) prefer to name "freezing index" the variable w 
defined by (II. 6p . having also in mind applications to frost propagation in porous media. 

Another meaningful variable of the Stefan problem is the enthalpy e, which in our case 
is defined by 

e = —dg"^ (— the variational derivative of \1/ with respect to 9), 






if u = -l 





if n < 1 





if u = +1 
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whence e = 6 + u = wt + u. Then, the governing balance and phase equations are given by 

et + divq=/ (1.8) 

ut + du'^ = (1.9) 

where q denotes the thermal flux vector and (i„^ stands for the variational derivative 
of \1/ with respect to u. Hence, (II. 9p reduces exactly to (II. 2p along with the Neumann 
homogeneous boundary condition for u. If we assume the classical Fourier law q = — V^ 
(for the moment let us take the heat conductivity coefficient just equal to 1), then (II. 8p is 
nothing but the usual energy balance equation as in the Caginalp model [5]. This is also 
as in the weak formulation of the Stefan problem, in which the mere pointwise inclusion 
u E (9/[_i_+i]) (9), or equivalently 6 E 9/[_i^+i](m), replaces (II. 2p . 

Another approach, which is by now well established, consists in adopting the so-called 
Cattaneo-Maxwell law (see, e.g., [U [H] and references therein): such a law reads 

q + £:qt = -V6', for 5 > small, (1.10) 

and leads to the following equation 

eOtt + e- Ae + eutt + Ut = f in ft x {0, T) 

which has been investigated in [H]. On the other hand, if we solve (ll.lOp with respect to 
q we find 

q = qo + A; * V6', where {k * V6')(x, t) := / k{t - s)V9{x, s)ds, 

Jo 

qo{x,t) is known and can be incorporated in the source term, k{t) is a given kernel 
(depending on e of course): from (II. Sp we obtain the balance equation for the standard 
phase fleld model with memory which has a hyperbolic character and has been extensively 
studied in [U [S] . 

In [HI [HI [ini [H] Green and Naghdi presented an alternative approach based on a 
thermomechanical theory of deformable media. This theory takes advantage of an entropy 
balance rather than the usual entropy inequality. If we restrict our attention to the heat 
conduction, these authors proposed three different theories, labeled as type I, type II and 
type III, respectively. In particular, when type I is linearized, we recover the classical 
theory based on the Fourier law 

q = -aVwt, a > (type I). (l-H) 

Furthermore, linearized versions of the two other theories yield 

q = -(3Vw, (3>0 (type II) (1.12) 

and 

q = -aVwt - (3Vw (type III). (1.13) 

Note that here we have used the thermal displacement (II. 6p (instead of 6) to write such 
laws. We also point out that (11.120 - 01.130 have been recently discussed, applied and 
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compared by Miranville and Quintanilla in [151 [IHl E] (there the reader can find a rich 
hst of references as well). In particular, fll.l3p leads via fll.Sp to our equation f ll.ip : further, 
a no flux boundary condition for q corresponds to dnW = in (11.31) . 

Thus, the system ([III])-([riD results from ([TBD-dlSD when ([13]) and flTT^ are postu- 
lated. We are interested in the study of existence, uniqueness, regularity of the solution to 
the initial-boundary value problem f ll.ip -( lL^ when 7 is an arbitrary maximal monotone 
graph, possibly multivalued, singular and with bounded domain. Of course, the case of 
\1/ shaped by a multiwell potential u 1— )■ —Wtu + 0(m) + G{u) is recovered as a sample. 
Then we study the asymptotic behaviour of the problem as /3 \ 0, obtaining convergence 
of solutions to the problem with (3 = 0, which corresponds to fll.lip . the (type I) case 
of Green and Naghdi. We also prove two error estimates of the difference of solutions in 
suitable norms, showing a linear rate of convergence in both estimates. In a subsequent 
study we would like to address the investigation of the analogous limit a \ to obtain 
the (type II) case in (11121) . 

The paper is organized as follows. In Section [2] we state the main results related to 
the problem fll.ip -f lL^ : existence and uniqueness of a weak solution, regularity results 
yielding a strong solution, further regularity results ensuring the boundedness oi u, Wt 
and of the appropriate selection oi'~f{u). Section [3] contains the related statements. Then 
we investigate the asymptotic limit as /3 \ 0: precisely, the convergence result and the 
error estimates under different assumptions on the data. In Section |4| we introduce some 
notation and present the uniqueness proof. The approximation of the problem (11.10 - 
(II. 4j) via a Faedo-Galerkin scheme and the derivation of the uniform a priori estimates 
are carried out in Section Regularity and boundedness properties for the solutions 
are proved in Sections [6H8l Finally, the details of the asymptotic analysis as /3 \ are 
developed in Section [91 



2 Well-posedness and regularity for a, /? > 

We point out the assumptions on the data and state clearly the formulation of the problem 
and the main results we achieve. Let fi C M^ be a bounded smooth domain with boundary 
T = dn and let T > 0. Set Q := fi x (0, T). We assume that 

a, /3 e (0,+cx)) (2.1 

/ e L'(0, T; H\ny) + L\0, T; L^Q)) (2.2 

7 C M X M is a maximal monotone graph, with 7(0) 3 (2.3 

(j) : M — y [0, +00] is convex and lower- semicontinuous (2.4 

0(0) = and 90 = 7 (2.5 

g : M — > M is Lipschitz-continuous (2.6 

Wo G H\n) , vo G L\n) , no G L^{n) , 0(no) G L\n). (2.7 

The effective domain of 7 will be denoted by Dij). We consider 



(2.12) 
(2.13) 
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Problem (P^^^). Find [w, u, ^) satisfying 

w e W^'°°{0, T- L'^in)) n H\Q, T; H\n)) (2.8) 

WtteL\0,T;H\ny) (2.9) 

u e H\0, T; H\Q)') n C° ([0, T]; L\Q)) n ^^(0, T; H\Q)) (2.10) 

^ G L2(g) , u E D{-i) and ^ G 7(m) a.e. in Q (2.11) 

for all V G if^(l]) and a.a. t G (0,T) 
(ui(t), i;) + (Vu(t), VtOi2(^) + (^(t), w)i2(f^) + (^(^)(t), t')L2(f^) = {wtit), tOi2(n) 

for all w G i/^(fi) and a.a. t G (0,T) 
w(0) = Wo in /7^(fi) , Wi(0) =t;o in i^^(^)', m(0) = Mq in i:^(fi). (2.14) 

We can prove the well-posedness of this problem. 

Theorem 2.1 (Existence and uniqueness) . Let assumptions (I2.ip - (l2.7p hold. Then 
Problem {Pa,^) has a unique solution. 

Next, in addition to (I2.ip -( ET|) . we suppose 

/ G L2(0, T; L^{n)) + L^(0, T; if ^(1^)) (2.15) 

Wo G ii2(fi) , dnWo = on r , VqE H\n) , Uq G if^(fi) ; (2.16) 

in this case, we are able to prove a regularity result, which allows us to solve a strong 
formulation of Problem (Pq,,/?). 

Theorem 2.2 (Regularity and strong solution). Assume (I2.15p - (l2.16p in addition 
to (I2.ip - (l2.7p . Then the unique solution [w, u, ^) of Problem {Pa,p) fulfills 

w G iy^'~(0, T; i/^(l^)) n H\ld, T- H\n)) (2.17) 

Wtt GL^(0,T;L2(fi)) (2.18) 

u G i/^(0, T; L\Q)) n C° ([0, T]; H\Q)) n ^^(0, T; ii^(fi)) . (2.19) 

In particular, {w, u, ^) solves Problem {Pa,i3) 'if^ « strong sense, that is, w and u 
satisfy 

Wtt — aAwt — f3Aw + Ut = f a.e. in Q 

Ut — An + C, + g{u) = Wt, { G 7(n) a.e. in Q 

dnW = dnU = a.e. onT X (0,T) . 

The aim of the subsequent results is to provide L°° estimates. We will need to 
strengthen again the hypotheses on the initial data. For s G D{'~f) let us denote by 
7^(5) the element of 7(5) having minimal modulus. Then, we require that 

Mo e H'^{Q) , dnUo = on r (2.20) 

Mo e D(7) a.e. in n , 7°(mo) e L'^in) . (2.21) 
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Theorem 2.3 (Further regularity). If the conditions ( ^2A} - (K7^ . ( l235|l - (12776]) and 
fl2.20p -f l2.2ip hold, then the solution [w, u, ^) of Problem {Pa,^) fulfills 

u G iy^'~(0, T- L\n)) n H\0, T; H\n)) n L°°(0, T; H^n)) . (2.22) 

The above results still hold if the dimension A^ of the domain ^2 is arbitary. On the 
other hand, since f l2.22p implies in particular that u is continuous from [0, T] to the space 
H^[Q) for all s < 2, then, if we let A^ < 3 and s sufficiently large, it turns out that 
H^{Q.) C C°(r2) and consequently 

u e C\Q) . 

Finally, we assume for the data enough regularity to get L°° estimates for Wt and ^. 
The hypothesis A^ < 3 is essential in the proof of the following result. 

Theorem 2.4 (L°° estimate for Wt and ^). In addition to assumptions f l2.ip - fl2.7p . 
fl27[5|) -f l216|) and (^^-^Mj, we ask 

f G L°°(0, T; L^in)) + ^"(0, T; H^{n)) for some r > 4/3 (2.23) 

7°(no)GL-(fi). (2.24) 

Then we have 

wtEL'^iQ), eeL-(g). 

Remark 2.5. All the statements contained in this paper still hold if fi C M'^ is, for instance, 
a convex polyhedron, for which standard results on Sobolev embeddings and regularity 
for elliptic problems apply. 



3 Asymptotic behaviour as /3 \ 

Let us fix the parameter a once and for all. We shall concentrate on the asymptotic 
behaviour of the solution as /3 \ 0, so we let /3 vary in a bounded subset of (0, +00). We 
allow the source term and the initial data in Problem (Pq-,^) to vary with /3, by replacing 
/, Wo, vq and uq in fl2.12p and fl2.14p with /^, wo,/35 'yo,/3 and mo,/3 respectively. We will 
denote by (w/3, m^, ^^) the solution to Problem (Pa,/?). 

If we set /3 = in the statement of Problem (Pq,,/?), we get a first-order system 
of differential equations, with respect to time, in the variable Wt, which is of physical 
relevance (recall that Wt = 9). Anyway, we avoid this change of variable, in order to 
preserve the formalism. We introduce the formulation of Problem (P^), in which /3 is set 
to be zero. 

Problem (Pq). Find {w, u, ^) satisfying fl2.8|) - fl2.11|) as well as 

{wu{t), v) + a (Vwtit), Vv)^,^^^ + {ut{t), v) = {fit), v) 
for aU V G H\n) and a.a. t G (0,T) 
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{ut{t), v) + (Vn(t), Vi;)i2(n) + ((^ + 9{u)){t), tOi2(n) = (w^iW, ^Ol^cq) . . 

for all veH\n) and a.a. te(0,r) 

w(0) =m;o in i/^(fi), Wi(0) = t;o in i/^(fi)' , m(0) = mq in ^^(fi). (3.3) 

We state at first the well-posedness of Problem (P^) and a convergence result. 



Theorem 3.1 (Well-posedness for (P^)). If the hypotheses (I2l2|) - (12171) hold, then 
Problem (Pa) admits exactly one solution. 



Theorem 3.2 (Convergence as /3 \ 0). We assume 02.2p -f l?7r|) and 

fp^f mL\0,T;H\ny) + L\0,T;L\n)) (3.4) 

Wo,/3 -^ Wq in H^{Q) , vq^j3 -^ vq , uq,/? ^ uq in L^(f2). (3.5) 

Then, the convergences 

wp^*w zniy^'~(0,T;L2(fi)), wp ^ w m H\0,T; H\n)) 

up^u m H\0, T; H\n)') n ^^(0, T; H\n)) 

/ioW, where {w,u,C,) denotes the solution to Problem (Pa)- 

With slightly strengthened hypotheses, we are able to prove the strong convergence 
for the solution and even to give an estimate for the convergence rate. 

Theorem 3.3 (First error estimate). In addition to 02.3p -( ETB|l and 03.4p - 03.5p . we 
assume 

\\f/3 - /llL2(o,r;Hi{f^)')+-f'MO,T;L2(n)) - '^ ^ (3-6) 

|ko,/3 - Woll^i(f^) + |K^0,/3 - Vo\\Hi(Qy + ll"0,/3 " Mo|li2(f^) < c /3 (3.7) 

for some constant c which is independent of (3. Then one has the estimate 

11^/3 ~ '"^ll/fi(0,T;L2(n))nL°°(0,T;Hi(n)) , . 

+ IF/3 — '"llL°o(0,T;L2(n))nL2(o,T;/fi(f^)) — ^^ 

where c does not depend on (5. 

If 7 is a (single-valued) smooth function, and if enough regularity on the data is 
assumed, it is possible to obtain much stronger estimates. The assumption A^ < 3 on the 
spatial dimension is essential for the proof of the following result. 
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Theorem 3.4 (Second error estimate). Let ^^-^M), (ED-lESD hold and 

7 : D{^) — > M be single-valued and locally Lip schitz- continuous. (3-9) 

Moreover, assume that the data {fj^, Wq,/?, "^o,/?! ""o./^}; O'S well as {/, Wq, Vq, Uq], 
satisfy f l2J[5|) - fl2:T6|) . ([22QD -(jMID , (E23])-(l221 along with 

ll//3|lL°°(0,T;L2(f7))+L'-{0,T;Hi(f^)) "^ \\'^^A\m{n) + \\1K^Q,I3)\\ L'^{p,) < C (3.10) 

II//3 - /llL2(o,r;L2(n))+Ll(0,T;Hi(!^)) - '^ ^ (^-H) 

||wo,/3 - Wo||^2(f^) + \\vo,i5 - ^oHhh^) + 11^0'^ ~ "o|li:^i(n) < c/3 (3.12) 

where r > 4/3. T/ien the estimate 

+ \\Ui3 ~'"ll//i(o,T;L2(Q))nL°o(0,r;Hi(n))nL2(0,T;//2(f7)) < ^p 

holds for a suitable constant c, which may depend on a but not on /3. 



(3.13) 



4 Notation and uniqueness proof 

Before facing the proof of all the results, for the sake of convenience we fix some notation: 

Qt = nx {0,t) for < t < T, Q = Qt, 

H = L^{n), V = H\n), W ={veH\n): dnv = a.e. onr}. 
We embed H in V', by means of the formula 

{y, v) = {y, v)^ ioT aWy e H , V e V . 

Furthermore, the same symbol ||-||^ will denote both the norm in L'^{Q) and in L'^{Q)^; we 
behave similarly with ||-||y^. If a, b are functions of space and time variables, we introduce 
the convolution product with respect to time 

{a*b){t)= a{s)b{t-s)ds, 0<t<T. 
Jo 

We also point out that the symbols c, c, - even in the same formula - stand for different 
constants, depending on fi, T and the data, but not on the parameters a, (3. However, 
as we will be interested in the study of convergence as /3 \ 0, if a constant c depends on 
a, /3 in such a way that c is bounded whenever a, /3 lie bounded, then we will accept the 
notation c. A constant depending on the data and on a, but not on /3, may be denoted 
by Cq or Cq, j or simply c, as it will happen in Section [91 

In our computations, we will often exploit the Holder and Young inequalities to infer 
/ ab<^ f \\a{s)\\lds+'^ f \\b{s)\\lds 

jQt ■^'^ JO ^ Jo 
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where a, b & L'^iQ) and a > is arbitrary. We point out another inequahty which will 
turn out to be useful: if v? G H^{0,T] H), then the fundamental theorem of calculus and 
the Holder inequality entail 



(/.(0)+ / ^t{s)ds <2||v.(0)||^ + 2T / yt{s)\\lds (4.1) 

^0 H Jo 



for all < t < T. Now, let us concentrate on the uniqueness proof. 

Let {wi, Ui, ^i) and {w2, U2, (,2) be solutions to the Problem (Pa,/3); we claim that 
they coincide. Setting w = Wi — W2, u = ui — U2 and ,^ = ,^1 — ,^25 we easily get 

{wtt{t), v) + a {Vwt{t), Vv)h + 13 (Vu'(t), Vv)^ + {utit), v) = (4.2) 

{ut{t), v) + (VM(t), V^)^ + (e(t), v)h + {g{ui){t) - g{u2){t), v)^ = {wt{t), v)^, (4.3) 
for all V E V and a. a. < t < T, along with the initial conditions 

w{0) = WtiO) = u{0) = . (4.4) 

We choose v = u{t) in equation 04. 3|) and integrate over (0,t); thus, we obtain 



\\Ht)\\l+ f\\Vu{s)\\lds+ I iu 

^ JQ JQt 



{g{ui) - g{u2))u+ / wtu. 



Accounting for the Lipschitz-continuity of g, the Holder inequality and the monotonicity 
of 7, frow the above equality we easily derive 

l\Ht)\\l+ [ \\\/uis)\\lds<c [ \\uis)\\lds+ I WtU. (4.5) 

•^ ^0 Jo JQt 

Integrating in time the equation (14. 2 p (this is possible thanks to (12.81) ) and taking the 
initial data (14.41) into account, we have 

{wt(t), v)^ + a iVw{t), Vv)j, + /3 (1 * Vw{t), Vv)^, + (^(t), t;)^ = ; (4.6) 

we choose v = Wt{t) in (14.61) and integrate over (0,t). Noticing that the equality 

(1 * Vwit), Vwtit))^ = Jt(^* ^^W' ^^W)h - II V^(t)|li (4.7) 

holds, we get 

ft 

J^ \\wt{s)\\lds + I \\Vw{t)\\l = -/3 (1 * Vw{t), Vw{t))^ 



+/3 / \\'Vw{s)\\fjds — / uwt^ 

Jo JQt 



(4. 



The Holder inequality and (14. ip allow us to deal with the right-hand side of this formula: 

?2 rt 
'0 



/3(l*Vw7(t), Vw7(t))^<^ / \\Vw{s)\\lds+'^\\Vw{t)\\l ■ (4.9) 

a Jo 4 
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Collecting now (gSD, (BSD and (gj]), it follows that 

\ \Ht)\\l + [ \\'^n{s)\\lds + f^ \\w,{s)\\lds + ^ ||V^(t)||^ 

<c I \\u{s)fHds + c((3 + —\ I ||Vw(s)||5^rfs; 

then, by applying the Gronwall lemma and recalling (I4.4p . we obtain u = w = Q almost 
everywhere in Q. A comparison in (12.131) and the density of H^{Q) as a subspace of L'^iQ) 
entail ^ = almost everywhere in Q; thus, the proof of uniqueness is complete. 

5 Approximation and a priori estimates 

We are going to prove the existence of a solution to Problem (Pq,^) via a Faedo-Galerkin 
method. First, we approximate the graph 7 with its Yosida regularization: for all e G (0, 1] 
say, we let 



le 



-{/ — (/ + £7) } and 4>e{s) := min S — |t — s| + 0(r) > for s G 



where I denotes the identity on M. We recall that 0e is a nonnegative, convex and 
differentiable function, 7^ is Lipschitz-continuous, monotone and 

7.(0) = 0, 0^ = 7., O<0,(s)<0(s), \le{s)\<\l°{s)\ V£>0, sGM (5.1) 

(see, e.g., [21 Prop. 2.6, p. 28 and Prop. 2.11, p.39] or P pp. 57-58]). 

We look for a solution of the approximating problem in a finite-dimensional subspace 
Vn '^ V, chosing a sequence {Vn} filling up V; then we get a priori estimates and use 
compactness arguments to take the limit as n — > +00. In a second step we let e \ 0. 

A special choice of the approximating subspaces will be useful. Let {fjjjgf^ be an 
orthonormal basis for V satisfing 

—Avi = XiVi in Q, dnVi = on F (5-2) 

where {Ajj^gj^ are the eigenvalues of the Laplace operator; also, let Vn be the subspace of 
V spanned hy vi, . . . , Vn, for all n &N. Thus, we have defined an increasing sequence of 
subspaces, whose union is dense in V, and hence in H; furthermore, we notice that the 
regularity of Q implies Vn C W, for all n G N. 

As approximations of the data wq, vq, uq we choose the projections on Vn- let wo,n 
be the projection of wq, with respect to V, and let wo,n, uo,n be the projections of vq, uq, 
with respect to H. We notice that 

Wo,n — > Wo inV , VQ^n — > vo in H , %,„ — > uq in H . (5.3) 

We also need to regularize the source term /: so, we first write 

/ = /(I) + /(2) , where f'^ G ^^(0, T; V) and f^^ G ^^(0, T; H) , (5.4) 
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then we assume /n , /« to be functions in C° ([0,T]; V), C° ([0,T]; H) respectively, 
such that 

/(i)^/(i) inL2(0,r;O, /f -^/('MnL^(0,T;i7); (5.5) 

we also set /„ = fi^' + fir . 

Now we are ready to state the approximated problem. For the sake of simplicity, we 
do not specify explicitly the dependency on e in the solution. 

Problem (Pq,,^)„ ^. Find T„ G (0,T] and {wn.Un) satisfying 

Wn G C2([0,T„]; Vn) , n„ G Ci([0,T„,]; K) 

(a>„(t), t;)^ + a (V9iM;„(t), S/v)^ + /? (VM;„(t), Vt;)^ + {dtUn{t), v)^ 

= {fn(t), v) for all V G Vn and all t G [0,T„] 

{dtUnit), v)jj + (Vn„(t), Vf)^ + (7e(u„)(t), t^j^^ + (5f(n„)(t), i;)^ 
= (dtWnit), v)jj for all f G K. and all t G [0,T„] 

W„(0) = Wo,n , dtWni^) = Vo^n , Mn(0) = Mo,n • (5.8) 

Writing Wn and u„, as linear combinations of fi, . . . , w„ with time-dependent coefficients, 
and testing equations (15. 6p and (15. 7p by t' = fi, . . . , w„, we obtain a system of ordinary 
differential equations, for whose local existence and uniqueness standard results apply. 
Thus, Problem (P^,/?)^ admits a solution, defined on some interval [0, T„]. The following 
estimates imply that these solutions can be extended over the whole interval [0, T]. 

First a priori estimate. We choose v = Un{t) in equation (15. 7p and integrate 
over (0, t): 



(5.6) 

(5.7) 



2 



Un{t)\\H+ / \\VUn{s)\fjjds+ / 'ye{Un)Un 
Jo JQt 

= - / g{Un)Un+ / UndtWn + - \\uo^n\\% ■ 
JQt JQt ^ 

The last term in the left-hand side is non negative, because 7^ is increasing and 7£(0) = 0; 
it will be ignored in the following estimates. Meanwhile, the right-hand side can be easily 
estimated using the Lipschitz-continuity of g and (15. 3p : so we get 

o ll^n(^)llH+ / \\^Un{s)\\\ds<c j \\un{s)\\\ds+ I UndtWn + C. (5.9) 

^ Jo Jo JQt 

Following the same computation as in the uniqueness proof, we integrate equation (15. 6p 
with respect to time: 

(dtWnit), v)h + a (Vwnit), Vv)h + /3 (1 * Vwn{t), Vv)^ + (M„(t), w)^ 

= (1 * fi^\t), V)+{1* f^\t), V)^ + iVo,n + Uo,n, v) ^ + « (Vu^Q.n, Vv)^ 

for all V & Vn and < t < T„. We take v = dtWnit) in the previous equation and integrate 
over (0,t). Recalling the identity (14.71) . we have 

/t 7 

\\d,w^{s)\\lds + I WVwnml = J2T,{t) + I ||V^o,„,|li (5.11) 
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where we have set 

Tl{t)=l3 [ \\VWn{s)\\lds, T2{t) = -(3{l*VWn{t),VWn{t))H 

Jo 

Ts{t) = - fundtWn, mt)= f{l*fi'\s),dtWn{s))ds, n{t) = / (l * /f ) 9*^, 

jQt Jo JQt 

Teit) = {vo,n + uo,n, dtWn{s))jjds, Tj{t) = a {Vwo^n,'^dtWn{s))jjds. 
Jo Jo 

We do not need any estimate on terms Ti and T3. With simple apphcations of the Holder 

inequality, we estimate T2, T5 and Tq: 



T2(t)<^||Vw„(t)||?, + — f \\\/Wnis)\\lds 

m)<\I^Hw^{s)\\lds + lj\l*fi'\s)\\lds 

Teit) < - / \\dtWn{s)\\%ds + c\\Vo^n\\% + c\\uo^n\\H ■ 

4 Jo 
We deal with T7 by direct integration and the use of the Holder inequality: 

o 

Now we pay attention to T4 and integrate by parts in time: 

Ut) = (1 * fi'\t), w4t)) - f {fi'\s), w4s)) rf. < ^ II 1 * fi'\t)\\l, 

Jo '^^ 

where a > is arbitrary, to be set later. According to the definition of the norm in V 
and the inequality (14. ip . we have 

1 ft 

Ut)< — \\l*fi'\t)\\l, + aTj^ \\d,w^(s)\\lds + ^\\Ww^it)\\l 



t rt / rs 



^ f \\f^'Hs)f 

I \\Jn \^J\\v^ 

'0 JO \J0 



+ 2 I \\fi'\s)\\^,ds + T I (I \\d,w^{T)\\l dr ] ds 



1 /■* 

-J \\Vwn{s)\\lds + T{a + l)\\wo,n\ 



2 

H ■ 



We collect all the terms containing ||t?tifn|lL2(oiH) ^^^ ll^^n(^)||// in the left-hand side 
of (15. lip : their coefficients turn out to be, respectively, 

ki = Ta , k2 = - i a] . 

2 ' 2 V2 / 

We choose a < min {a/4, 1/4T}, so that ki > 1/4, k2 > a/8. We also remark that the 
assumptions (15. 5 p and (15. 3 p enable us to get a bound for terms involving /„ , fn and 
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the initial data. Finally, adding fl5.9p and (15. lip and taking into account all the previous 
inequalities, we obtain 

1 2/ 2I/ 2 (^ 1 

\ i \\dtWn{T)\\H dr ] ds + Ca \\VWn{s)\\H ds + Ca . 



< c I \\uJs)\\%ds + T 



t l-t 

( "lll^ 

n\^)\\H 
'0 

The Gronwall lemma entails 



Second a priori estimate. Since 0^ is at most of quadratic growth, by definition, 
and 7£ is Lipschitz-continuous, from the estimate fl5.12p we directly derive 

l|0eK)||ico(o,T;Li(f7)) < cL,l (5.13) 

ll7eK)||i2(Q) <42; (5-14) 

where the symbols c^ ^ denote positive constants, possibly depending on e and a, but not 
on n and /3. 

By ()5.2p . we can easily check that 

iy^ z)h = (PriV, z)h for all y G l^ , 26 K 

where PnH is the projection of y in V^, with respect to V . Then, as we have a uniform 
estimate for m„ in L^(0,T; V), it is not difficult to extract from fl5.7p the property 

||5iMn||L2(o,T;y') -^°.3- (5-15) 

Third a priori estimate. We take v = dtWn{t) as a test function in equation fl5.6p 
and integrate over (0,t); thanks to the Holder inequality, we get 

I Wdtw^ml + « ^ \\VdtWn{s)\\l ds + ^ ||V^„(t)||?, 

< f\fi'^-dtUnis),d,Wnis))ds+ f\\fi'\s)\\^\\dtWr,,{s)\\^ds (5-16) 

Jo Jo 






We consider the term involving /„ — dtUn'- 



/ {fn^ -dtUnis), dtWnis))ds<- \\fn\s)\f ds + - \\dtUnis)fy, ds 

Jo (^ Jo «io 

2 a [^ 2 

^ Jo 



a '" 



Because of the estimate f lS.lSp and the properties fl5.5p and fl5.3p . from fl5.16p we deduce 

■^ Jo Jo 
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where d depends on e, a. Hence, by a generalized version of the Gronwall lemma (see, 
e.g., [21 pp. 156-157]), we infer that 

||^n|liyi,oo(0,T;_f/) ^ y^ ll^"ll//i{0,T;V) — '^a,4 • l^-l' J 

Passage to the limit as n — )■ +oo. From the estimates ( 15.1 2p . ( l5.13p -( 1B.15p . ( I5.17p . 
with standard arguments of weak or weak* compactness we can find functions (w^, u^ 
such that, possibly taking a subsequence as n — > +oo. 



t'n - 


^*w. 


in l^i'°°(0,T;i7)nL°°(0,T;^) 


(5.18) 


Wn - 


-^ We 


in H^{Q,T-V) 


(5.19) 


Un 


-^Ue 


in H\Q,T-V')fM^{id,T-V) 


(5.20) 


Un - 


-*Ue 


in L°^{0,T-H). 


(5.21) 



Note that (15.190 implies the strong convergence 

w^^w, inC°([0,T];if); (5.22) 

on the other hand, the generalised Ascoli theorem and the Aubin-Lions lemma (see, e.g., 
PSI pp. 57-58] and [IHl Sect. 8, Cor. 4]) entail 

n„ — > u, strongly in C^ ([0,T]; V') and in L^{Q)] (5.23) 

thus, since g and 7^ are Lipschitz-continuous, we easily check that 

giun) — > giue) and ^leiun) — > L strongly in L^(Q), 

where ^^ = 7e(Me). We then take the limit as n — > +00 in (I5.6p - (l5.8p and see that 
(wg. Me, ^e) fulfills equations (12.110 - 02.140 . where 7 is replaced by 7^. Indeed, by (15.220 - 
(I5.23P and (15. 3p . it is obvious that We(0) = Wq, u^iQ) = Uq. To deal with the last initial 
condition properly, we fix a test function v G Vm-, where ?n > 1 is arbitrary, and we 
integrate in time equation (15.60 : we get equation (I5.1O0 . for < t < T and n > m. 
Arguing as in [131 PP- 12-13], we can take the limit in (I5.1O0 . (15.70 and check that 
(w^£, ^e, is) fulfills 

{dtw.it), v) = -a iyw.it), Wv)h - /3 (1 * Ww,{t), Wv)^ 
- (ueit), v) + {I* fit), v) + a (Vwo, Wv)^ + {vq + Mo, v)^ 

(dMt), v) + (Vn,(t), Vm)^ + iUt), v)j, + {g{u,){t), v)j, = {dtw,{t), v)j, (5.25) 

for a.a. t G (0,T), m > 1 and v G Vm', by a density argument, the same equalities hold 
when V &V. Since the right-hand side in (I5.24p is a continuous function in [0,T], taking 
t = we find that 

{dtWs{0), v) = [vo, v)^ for all m G ^ 

whence the second of (I2.14p follows. 

Fifth a priori estimate. As a consequence of the weak lower semi-continuity of the 
norm in a Banach space, {we, Ue, C,e) satisfy the estimate (I5.12p : we now need to improve 
estimates fl57[3D - (15351) . (IHTTjl . 
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We first notice that, because of tlie Lipscliitz-continuity of 7^, C,e{t) G V for all t; thus, 
we can choose v = ^^(t) in equation fl5.25|] and integrate over (0,t), to get 

/ dtu,^,+ [ j',{u,)\Vu,f+ [ \\^,{s)\\lds= [ g{u,)^,+ [ dtw.i,. (5.26) 
JQt JQt Jo JQt JQt 

In view of fl5.ip . we have 

f 9 

on the other hand, because of the Lipschitz continuity of g, 

9iUe)^e<cf (K| + l)ee<c/ {\\uM\\l + '^)ds + l f \\Us)fHds. 

From these estimates and fl5.26p . we derive 

I Unem+ I i,{u,)\Vu^' + \ f\\Us)\\\ds 
Jn JQt ^ Jo 

<c \\ue{s)\\lj ds + c / \\dtWe{s)\\lj ds + / (f)s{uo) + C. 
Jo Jo Jn 

We notice that the second term in the lef-hand side is nonnegative, because of the mono- 
tonicity of %. Secondly, accounting for fl5.12p . (15. ip and (12.71) . we infer that 

\\MUe)\\Lo^(0,T;LHn)l) + he{Ue)\\L2(Q) < Ca ■ (5.27) 

Now, by comparison in the equation (I5.25p . we have 

ll^t'"e|lL2(o_T;V^') — '^a ', (5.28) 

and consequently we can also establish the estimate (I5.17p . now for a constant which is 
independent of e. 

Passage to the limit as £ \ 0. We are able to repeat the compactness argument 
as above and find {w, m, ^), a candidate for the solution to Problem (Pa,^), as a limit of 
a subsequence of {we, u^, ^e). The proof will be easily completed by the passage to the 
limit as e \ 0, provided that we deduce (12. lip . 

By construction, we can assume that 
from which the equality 



^\o Jq Jq 

follows; at this point, we apply P, Prop. 1.1, p. 42] and deduce (12. lip . Thus, the proof of 
the existence of a solution to Problem (Pa,/?) is complete. 
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6 Regularity and strong solutions 

This section is devoted to the derivation of further a priori estimates on the approximating 
solutions {wn, Un, ^n), which are independent of n and e, under stronger assumptions. 
The same compactness - passage to the hmit arguments then apply, and this will prove 
Theorem 12.21 We first notice that the hypothesis fl2.16p and Vn^W make it possible to 
assume 

"W^o.n — > Wq in VT, t>o,„, — )■ Vq and %,« — ^ Uq 'm.V ] (6.1) 

on the other hand, owing to fl2.15p . we can require fn G L'^{Q), fn € -^^(0, T; V) for all 
n G N and 

/(i)^/(i) inL^Q), /f ^/(^) mL\0,T;V). (6.2) 

Sixth a priori estimate. We choose v = dtWn{t) in the equation 05. 6p and integrate 
over (0,t); an application of the Holder inequality yields 



^ < - / dtUndtWn 



-\\dtWn{t)\\]j + a / \\WdtWn{s)\\\ds + -\\VWny^)\\H 

[ 1/3 

+ / \\fn{s)\\H\\dtWn{s)\\Hds + -\\Vo,n\\\ + -\\VWQ,n\\\ ■ 

Now, we take v = dtUnif) in (15. 7p and integrate over (0,t); recalling that 7^ = 0^, using 
the Holder inequality and the Lipschitz-continuity of (?, we get 

1 /"* 1 

7: / \\dtUn{s)\\\ds + -\\VUn{t)\\\+\\<i),{Un{t))\\^,,^. 
< / dtUn.dtWn + C (||n„(s)||^ + 1) ds+ ||0e(uo,n)||ii(Q) • 



Adding (16.31) and (16. 4p . thanks to the assumptions (12. 7p . (15. 3p . the inequality (14. ip and 
0£ < 0, we finally have 



-\\dtWn{t)fH + a / \\VdtWn{s)\\% ds + -\\Vwn{t)\\% 
\\dtunis)\\lds+l l|Vn„(t)||^ + WMMmimn) 



2 Jo 2 

<cj (J \\dtUn{T)\\ldT\ds + j \\Us)\\^\\dtWn{s)\\Hds + C. 

The generalised Gronwall lemma (see, e.g., 0, pp. 156-157]) enables us to achieve 

||'W^n||vi/i.«>(0,T;H) "I" V Q^ ll'^'"ll_ffi(0,T;y) ~^ V P ll*^'" llL°°(0,T;y) — ^^ V^-^) 

ll^"ll//i(0,T;H)nL°°(0,T;V) — ^2- (6.6) 

Remark 6.1. Only the hypotheses (I2.ip - (I2.7I) and / G L^{0,T; H) have been effectively 
exploited in the proof of this estimate. 
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Remark 6.2. By means of f l6.5p - fl6.6p . the estimates fl5.27p -f l5.28p can be rewritten in 
terms of some constant which is independent of a. 

Seventh a priori estimate. We take v = —Aun(t) in equation (15. 7p : this is possible, 
because of the special choice of the approximating space Vn- We integrate over (0,t) and 
use the Holder inequality and the Lipschitz continuity of g: 

- ||Vu„,(t)||^ + / \\Aun{s)\\jjds+ 7^(n„,) |Vu„| 

^ ^0 JQt 

= - g'{un) |Vn„|^ - / dtWnAUn + - || Vno.„||^ 
jQt JQt ^ 

< c||Vm„||^2(o,T;H) + 2 ll^*^"lli2(0,T;H) + g / \\^^n{s)\\H ds + - WS/Uq^uWh ■ 

The monotonicity of 7^ yields that the last term in the lef-hand side is non negative. 
Owing to conditions (16. ip on the data and estimates (16.51) . (16. 6p . we have 



1 /■* 

-/ \\Aun{s)\\\ ds < c forallO<i<T: 

2 Jo 



hence, on account of this inequality, the estimate (16. 6p and the boundary conditions for 
Uni known regularity results for elliptic problems entail 



\u 



"IIl2(0,T;VK) 



<C3 



(6.7) 



where C3 does not depend on a, /3. 



Eigth a priori estimate. Since Wn G C^([0,T]; Vn)-, the special choice of Vn enables 

us to take v = —AdtWnit) as a test function in the equation (15.61) . We integrate over (0, t) 
and use the Holder inequality: 

^ WVdtWnml + «/ \\mWnis)\\lds + | || A«7„(t) Hj, 
<^ f\\AdtWn{s)fHds + - f\\d,Un{s)\\lds + - t \\fi'\s)\\l ds (6.8) 

^ Jo "Jo "Jo 

- I fj,'^Ad,Wn + \ \\Vvo,n\?H + f II ^^^^0,^11^ ■ 



For the term involving /A, we integrate by parts in space, recalling that dnV = for all 

ve Vn. 



f^'^d, 



Wr, 



v/f ) ■ Vd, 



Wr. 



< f\\^fi'\s)\\^\\Vd,Wn{s)\\^ds 

Jo 



(6.9) 



Then, in view of (16.11) . (16. 2p . (16. 6p and owing to the generalized Gronwall lemma (see [2], 
pp. 156-157]), from ([S3)-(IS3D we obtain 



ll^nll^i, 00(0, T;y) ~^ V " ll'^"ll_f/i(0,T;VF) ~^ V P II'^"IIl°°(0,T;W) — '-'".4 • (6.1U) 
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Finally, if we choose v = dfwnit) in the equation (15.61) . we get 

thanks to the estimates above, it is easy to derive 

Having established all the a priori estimates corresponding to fl2.17l) - fl2.19p on the solu- 
tions of the approximating problem, we have completed the proof of Theorem 12.21 

7 Further regularity 

Throughout this section we assume (12.201) and (12.211) in addition to all the hypotheses 
we had in Section [61 As we are interested in proving Theorem 12.31 we should get further 
estimates on the solution of the approximated problem. By the stronger assumptions on 
the initial data, we can require 

Uo,n — >uo iriW. (7.1) 

Consider the equation (15. 7p and derive it, with respect to time, obtaining 

+ {g'{Un{t))dtUn{t), v)h = {d'fWnit), V)^ 

for all V EVn and a. a. t G (0,T). We choose v = dtUnif) as an admissible test function, 
integrate over (0,t) and use the Lipschitz continuity of g to get 

l^\\dtUn{t)fH+ I \\S/dtUn{s)\\\ds+ j i^{un)\dtuj^ <c j \\dtUn{s)\\\ds 

^ Jo JQt JO fY 2) 

r-t -j^ \ ■ ) 

+ / \\dl'Wn{s)\\^\\dtUn{s)\\^ds + -\\dtUn{Q)fH . 

Jo ^ 

Since the last term in the left-hand side is non negative because of the monotonicity 
of 7e, if we had a bound for the last term in the right-hand side, we could use the 
generalized Gronwall lemma to conclude. In order to provide such an estimate, we set 
t = 0, w = 9iM„(0) in the equation (15.71) : we obtain 

and thus, taking into account the Lipschitz continuity of (?, we infer 

||5t«ri(0)|| < ||Ano,„||^ + ||7^IIl-(m) ll%,n - ^o\h + II7sK)IIh 

+c(||no,„||^ + l) + llt^o.nlln- 

Now, assumptions (17. ip and (15.31) . as well as (12.210 and 17^1 < |7°|, enable us to achieve 

||5tM„(0)|| <c (7.3) 
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for all e > and n large enough, depending on e\ these requests on parameters are not 
restrictive, as we first take the limit for n — )■ +00, then for e \ 0. From fl7.2p and f l7.3p 
we deduce that 

||Mn|lvi/i.°°(0,T;_H')nii'i{0,T;V) — Ca,6 • ('-4) 

Finally, we consider equation f l5.7p and we rewrite it in the form 

(VM„(t), Vv)^ + (7eK(t)), v)h = {Fn{t), V) , 

for all V E Vn and a.a. t G (0,T), where F„ = dtWn — dtUn — g{un). Testing with 
V = —Aun{t) the previous equation and integrating by parts in space, we obtain 

l|A«n(t)||i + / 7^K(t)) \Vunit)f < ||F„(t)||^ ||An„(t)||^ for all < i < T . 
Jn 

Since the estimates (16.51) and (17.41) entail 

II F II < r 

ll-^n|lL°°(0,T;H) — ^ 

and we can apply the regularity results for elliptic problems, we deduce 



thus concluding the proof of Theorem 12.31 



8 L^ estimates 

The aim of this section is to obtain L°° estimates on Wt and on ^, under the hypotheses 
(K2^ and (Km . 

We first deal with Wt- Setting ip = awt + /3w, Theorem 12.21 entail that the equalities 

-ipt - A(p = -Wt -ut + f in Q , 5nV5 = on F x (0, T) 
a a 

hold almost everywhere. Furthermore, the assumption (I2.23p . the estimates (16. 5 p and 
(I7.4p and the continuous embedding V M- L^{Q) (valid if fi C M'^ is a bounded Lipschitz 
domain), yield 

-wt-ut + f e L°°(0, T- H) + L70, T; L^(fi)) , with r > 4/3 . 
a 

In these conditions. Theorem 7.1 in [T^ p. 181] applies and ensures that ^p G L°°{Q). 
Since we already know that w G L°°{Q) (as it is implied, for example, by (I6.10p ). we have 
Wt G L°°(Q) and 

1 c/3 

iRtllLoo(Q) < — IIV^IIloo(q) + — |r|Il°°(0,T;H/) — ^".8 • (8-1) 

We notice that, being a fixed and letting /3 vary in a bounded set, we can find an upper 
bound for the constant Cqs- 
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In order to prove a L°° estimate for ^, we consider the solution (w^, Ue) to the approx- 
imating problem, in which the Yosida regularization appears; we then fix p e (1, +00) 
and get a bound for ||7e(Me)||j;^p(Q), which is independent of p, e. From this, we will obtain 
a uniform bound for 

^^' p— >+oo ^^' 

and, via a weak* compactness argument, ^ G L°°{Q). For the sake of simplicity, we do 
not plug in the subscript e in the solution any more. 

We know that the equalities 

ut - Au + 'ys{u) + g{u) = wt in Q , (8.2) 

dnu = on r X (0, T) , u{0) = uq in (^ 

hold a.e.; we choose |7e(M)|^~ 1e{u) as a test function, by which we multiply both sides 
of the equation (18. 2p - this is admissible since u G L°°{Q). Integrating over Q, we get 

^0,,,(n)+ /vn-V(|7e(n)rS.(n))+ I \le{u)r' 



jQ 

where we have set 

0.,pW = / \le{s)r^le{s)ds for alH e M; 

JQ 

■js is increasing and 7e(0) = 0, so we have (f)e,p > for all e, p. Since Wt, u G L°°{Q) and 
g is continuous, for the right-hand side we have 



/ 

JQ 



'Wt-g{u))\-fe{u)f le{u) 



< Ca \\le{u) 



\Lp{n) 



on the other hand, a direct calculation and the monotonicity of 7^ show that 

VM-V(|7.(M)rSe(w)) =pi,{u)\^,{u)r'\Vu\^>0 a.e. in Q . 
Collecting all the information we have obtained so far, from (18.31) we derive 

'n ~ "' ~ " Jn 

and, since the first term can be ignored, we need only to find an estimate for the last 
term. We recall that, for the Yosida approximation of a maximal monotone graph, the 
inequality 

17^(5)1 < |7°(s)| for all s G D{-f) , e>0 

holds (see, e.g., [21 Prop. 2.6, p. 28]); according to that, we have 



(j)e,p{u{T)) + he{u)\\lp+i^Q) < Ca ||7£(«)IIlp(q) + / <i)e,p{uo) (8.4) 



n Jn P+^ Jn P+^ Jn 



P+^ Jn P 

P I 1 .0/.. nIp+i , -'- II.. iip+i 



< 



P 



I U/ XU'T'-^ I II llP+i 

T / b ("o)| +— Ylko|lLP+i(n) , 

t/ iZ IT ^ 
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where the Holder and Young inequahties have been used. We recall that uq G L°°{Q) by 
the assumption fl2.22p and also notice that the same inequalities imply 



Ca l|7eHlliP(Q) < — — l|7e(«)||ip+i(Q) 



p+1 '"'^ ^"^--nw p+i 

Now, we come back to the equation (18. 4p : according to the previous estimates, we infer 
that 

II / MIP+l ^ P II 0/ \I|P+1 I ■'■ II IIP+1 I '^a 

and, hence, 

ll7eHllLP+i(Q) < \P\n M|lLP+Hf^) + ll^o|lLP+i(n)+Ca| 

<Ca[ ||7°M Loo(n) + II%IIl-(q) + 1 
which provides the desired estimate and concludes the proof. 



9 Well-posedness of (P^^) and convergence as ^ \ 

Now we set the notation as in Section [3], since we are interested in the proof of Theo- 
rems I3.1H3.4[ We assume that the hypotheses (I2.ip - fl2.7p are satisfied, and we start by 
studying the convergence as /3 \ 0, by a compactness argument. 

Convergence as /3 \ 0. We recall the a priori estimates flHTTD . <KW\ . iKTf} . fl5:28|) 
which are independent of (3 and thus holding also for (w^, M/3, ^^). Moreover, adopting 
the notation as in fl5. 40 - 05. Sp . by a comparison in fl2.12p we find out that {d^Wi3 — f^ } 
is uniformly bounded in L^(0,T; V). Therefore, we can find a subsequence /3fc \ and 
functions w, u, ^ such that 

wp^^*w in W^''^{0,T;H), wp^^w in H^{Q,T-V) 

dlwp-ff^ - wu-f^^^ in L\Q,T-V') 

M/3j, tends to u weakly in i7^(0,T; V') fl L^(0,T; V^), whence strongly in L?{Q), 

ip, - i in L\Q) 

as k — ;■ +00, and here part of fl3.4p has been used. Then, in view of fl2.6p . fl3.4p and 
(13. Sp . we can pass to the limit in fl2.12p and (12.130 . as well as in the initial conditions 
(I2.14P which can be recovered weakly in V at least. On the other hand, u G -D(7) and 
^ G 7(m) a.e. in Q follow as a consequence of the above convergences and [H Lemma 1.3, 
p. 42]. 

Uniqueness for (Pq). By applying the previous result with /^ = /, WQ^p = Wq, 
vo,i3 = vq and mo,/3 = mq given, we obtain the existence of a solution to Problem (Pa); we 
still have to prove the uniqueness. Let {wi, ui, ^i) and {w2, U2, ^2) be solutions of (Pa); 
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we write down the equations for the differences w = wi — W2, u = ui — U2, ^ = ^i — $,2 
and integrate with respect to time the first one: 

{ut{t), v) + {Vu{t), Vv)h + (eW, v)h + (gi^iM - ^M(t), v)h = {wt{t), v)h , 

to be complemented with null initial conditions as in f l4.4p . We set v = Wt{t) in the first 
equation and v = u(t) in the second one, integrate over (0,t) and add the two equations; 
it is straightforward to obtain 

1^ \H{s)\\lds + I \\Vw,{t)\\l + \ Uml + [ \\^n{s)\\lds < cj^ \Hs)\\lds . 

According to the Gronwall lemma and owing to w{0) = 0, it turns out that w = u = 
a.e. in Q and, by comparison in the second equation, ^ = a.e. in Q. 

Error equations. Because of the uniqueness, the whole family {(if/^, up, C,i3)}a^Q 
converges, as /3 \ 0, to the solution {w, u, ^) of Problem (Pa)- So, it makes sense to 
study the speed of this convergence. In order to perform that, we set wp = wp — w, 
Ui3 = Up — u, ^13 = ^j3 — ^ and consider the problem obtained for these variables, by 
subtracting side by side the equations of Problems {Pa,p) and (Pa). For all v E V and 
a. a. t G (0,T), the equalities 



(d^w^it), v)+a iVdtw^it), Vv)j^ + (3 (Vw^(t), Vv)j^ + {dtUp{t), v) 



(9.1) 



(dtu^it), v) + (VM^(t), Vv)^ + Up{t), v) + (gMit) - g{u){t), v)^ 



H '-' "'^ ' ■-' '^ " '« (9.2) 

= {dtwp{t), v)^ 

are satisfied, as well as the initial conditions 

Wp{Q) = Wo,i3 , dtW(s{0) = Wo,/3 , M/3(0) = Mo,/3 , 

where fp = fp-f = jf^ + f^;\ 

Jf^=ff^-f'^ -^ inL2(0,T;O 
fi2) = fm_fi2) _^ i^ L\0,T-H) 

(cf. dM])), uJo,/3 ■= w^o,/3 - Wo, Wo,/3 := Vo,fi - vo, and Uo,/3 := tio,/3 - ^^o- 

First estimate for the convergence error. Now, we want to show Theorem 13.31 
so we assume all the needed hypotheses. Choose v = Uf^it) in the equation (19.21) and 
integrate over (0, t); by the monotonicity of 7 and the Lipschitz-continuity of g, we easily 
derive 

J- "^ '.^i|2 , / UT-,^ I Ml2 7 ^ ^ w^ Il2 , / 11^ / M|2 



^/3WL+ / l|Vn^(s)llHC?s< -||no,^||^ + c / \ufi{s)\^ds\ \ u^dtWf^. (9.3) 

^ Jo JQt 



GiACOMO Canevari and Pierluigi Colli 



23 



We integrate with respect to time the equation f l9.ip : 

{dtWf,{t), v)^ + a (Vw^it), Vv)h + /3 (1 * Vw^{t), Vv)^ + (up{t), v)^ 

= (1 * fl3{t),v) + {vq^p + Mo,/3, v)^ + a (Vwo,/3, Vw)^ . 

We set V = dtWp and integrate over (0, t); keeping only the first two terms in the left-hand 
side, we obtain 



\\dtwM\\\ds + - \\Vwp{t)\\\ < - WVwoAh 



(3{1* Vwi3{t), Vwfi{t))j^ + (3 (Vwfiis), Vwi3{s))j^ ds - up dtWp 



+ / {l*Jf\s)+ ^o,/3, dtwuisyj ds+ I (l*Jf^ + So,/3J dtwn + a VwQ^nVdtWp . 

(9.4) 



Qt 



Qt 



Due to the Young and Holder inequalities and the boundedness of {1^/3} in L^(0,T; V), 
we have that 



-13 (1 * vwp{t), vMt))^ < ^p' J^ \\vwp{s)\\l ds + ^ W'^Mml 



a 



(9.5) 



<cl3' + -\\VMmH 



and 



(3 I iVwp(s), Vidp(s))jjds < c(3^ + a \\Vwp(s)\\lds 
Jo Jo 

f ^22 

« / Vwo,pVdtW(3 < — \\Vwfi{t)\\^ + ca WVwo^/^Wfj . 
jQt -"-^ 



(9.6) 



(9.7) 



On the other hand, arguing as in the estimate of the term T^lt) of (15. lip we deduce that 

1 * fl^\^) + ^o,/3, dtWfi{s)j ds 
'l*jf\t)+vo,p,Mt))- / {fl'\s),wp{s))ds 



< c 



HI)/ 



fy{s) ^^ds+\\vo^p\\y, + \\wo^p\\HJ +- / \\dtwp{s)\\j^ds 



{9.1 



'0 
a ,, /.x,|2 , / / / II a ^ / Ml2 7 \ r , / WT-,^ I Ml2 



12 



l|Ve^(t)||^ + c/^(^^'||9,iI;^(r)||irfr^d5 + ca^||ViI;^(.)||^rfs 



Finally, we observe that 



l^f^^ ^u^^^diWf,<c[\f^^ 



L^{0,T;H) 



I^oMh 



\dtwp{s)\\lds (9.9) 
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Now we add fl9.3p and (19 ■4p : collecting also all the estimates in f l9.5p - fl9.9p . we find out 
that 



\ '\\dtWp{s)\\\ds+'^\\Vw,m\ 

^ ^0 ^ 



\\\'^M\l+ I \\^Ms)\\lds+^ 



<cl3' + c 



r(l) 
//3 



L'2{0,T;V') 



ft2) 
//3 



+ ll^0,/3||_f/ + ||^0,/3||v" + ll'"^0,/3||v^ 



+c / \\ui3{s)\\^ds + C 

'0 Jo \Jo 



L^{0,T;H) 

\dtWf3{T)\\j^dT \ ds + ca / || V'u;^(s)||^ (is. 

At this point, it suffices to recall (I3.6p - fl3.7l) and apply the Gronwall lemma to obtain the 
thesis of Theorem 13.31 



Second estimate for the convergence error. Our aim is to prove Theorem | 
whose hypotheses are assumed to be satisfied. Thus, we can apply Theorems 12.31 and 12.41 
to get a bound 



l'«/3lh 



\U\ 



\m, 



Uh 



< Cn 



(9.10) 



with Cq, which is independent of /3. Now, if 7 is a maximal monotone graph which reduces 
to a single- valued function in its domain, then .0(7) is an open interval (a, b) and, if 
b < +CXO, then 7(r) /^ +00 as r /^ 6; similarly, if a > — cxo then 7(r) \ —00 as r \ a. In 
any case, the condition (19.101) implies the existence of some compact interval K C ^(7) 
such that up{Q) C K for all /3 > 0, u{Q) C K. Since 7 is assumed to be locally Lipschitz- 
continuous (cf. (13. 9p ). thanks to (13.81) we immediately deduce that 



U^-iW 



L°°{0,T;H) 



< C ||M/3 — u\ 



L°°iO,T;H) 



<cl3. 



Moreover, by suitably modifying g we can set ^^3 = in equation (19. 2p . without loss of 
generality. 

We start by taking v = dtWjs in (19.11) . v = dtUp in (19. 2p . integrating both equations 
over (0, t) and adding side by side. Thanks to the Lipschitz-continuity of g and the Young 
and Holder inequalities, it is straightforward to obtain 

^\\dtwp{t)\\% + a J \\VdtW^{s)\\%ds + J \\dtUp{s)\\lds + -\\Vup{t)\\l 

<■ — lF/3|lL2{0,T;y) + T 

ft 



a 



. I ,\Vdtw^{s)\\%ds + c I \\ufs{s)\\%ds+- \\dtU/s{s)\\]jds 
^ Jo Jo ^ ^0 



+- 



a 



J/3 



L^{0,T;V') 4 



a 



L2{0,T;H) + 4 



\VdtWfj{s)\\j^ds 



I ff\s) ^\\dtw^{s)\\^ds + -\\vo,4% + -\\Vuo,p\\l . 



Taking into account conditions (13. 6p . (I3.12p and the previous estimate (13. 8p . we easily 
have 



\dtWl3 



2 a 
H+ 2 



\VdtW,s{s)f^ds 



ll\\dtMs)\\lds+^-\\VMml 



<cP' + 



ff\s) ^\\dtwp{s)\\^ds 



H 
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whence, by (13. 6 p and a generalised Gronwall lemma (cf., e.g., [21 Lemme A5, p. 157]), we 
infer that 

ll^/3lliyl.°°(0,T;ii')nHl(0,T;y) + \\'^P\\m{0,T;H)nL°°{0,T;V) — '^ P [^■^^) 

where the constant c obviously depends on a. 

Next, observe that the assumptions on the data are strong enough to guarantee that 
(19.11) and (19. 2p can be reformulated as 

d^wp — aAdtWjs = (3Awi3 — dtUp + fp a.e. in Q (9.12) 

dtUp - Aufs + g{up) - g{u) = dtWp a.e. in Q (9.13) 

along with the homogeneous Neumann boundary conditions for both Wj3 and up. 

In view of (19. lip , by a comparison of terms in (I9.13P it is standard to deduce that 
||^w^|l2,2(or-H) — C/9 and consequently, owing to elliptic regularity estimates, we obtain 

II^^IIl2(o,T;H/) ^Cc./3. (9.14) 



At this point, let us emphasize that for the proof of (19. lip and (19.140 we have just used 
the control (13.60 on the difference fp. 

We now pay attention to the equation (19.121) and multiply both sides by —AdtWj^, 
which belongs to L'^iQ) (cf. (I2.17P ). and integrate, also by parts, over Qt- By means of 
the Holder and Young inequalities, we infer that 



- \\VdtW,3 



H 



a 



/•* 2 1 2 /3^ /■* 2 

Jo ^ (^ Jo 



— / \\dtUfj{s)fjjds + - 
a Jo "a 



f(i) 



L^{0,T;H) 



+ 



a 



v/, 



7i^)i 



H 



\\AdtWi3{s)f^ds 
\VdtWfj{s)\\j^ds. 



Hence, recalling the uniform boundedness of {w/3} in L^(0,T; W), we use (I3.12p . (19.110 . 
(13. lip and apply the generalised Gronwall lemma as before to obtain 



WVdtW^m'^ + / \\AdtWi3{s)\\'^ds<cP^. 
Jo 

Now, by virtue of (14. ip and (I3.12p we also infer 

\\Awp{t)\\j^ < cf3 for all t e [0,r]. 



Then, standard elliptic regularity properties and the previous estimates (19.111) and (I9.14p 
lead us to (I3.13p . thus completing the proof of Theorem 13. 4[ 
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